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ABSTRACT 

Context. Evolution of galaxies through cosmic time has been widely studied at high redshift, but there are a few studies in this field at 
lower redshifts. However, low-redshifts studies will provide important clues to the evolution of galaxies, furnishing the required link 
between local and high-redshift universe. 

Aims. In this work we focus on the metallicity of the gas in spiral galaxies at low redshift looking for signs of chemical evolution. We 
analyze the metallicity contents of star forming galaxies of similar luminosities at different redshifts, we studied the metallicity of star 
forming galaxies from SDSS-DR5 (Sloan Digital Sky Survey-Data Release 5), using different redshift intervals from 0.1 to 0.4. 
Methods. We used the public data of SDSS-DR5 processed with the STARLIGHT spectral synthesis code, correcting the fluxes for 
dust extinction, estimating metallicities using the R 2 j method, and analyzing the samples with respect to the [N n] /16583/[0 n] /13727 
line ratio. 

Results. From a final sample of 207 galaxies, we find a decrement in 12+log(0/H) corresponding to the redshift interval 0.3 < z < 0.4 
of ~0. 1 dex with respect to the rest of the sample, which can be interpreted as evidence of the metallicity evolution in low-z galaxies. 
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1. Introduction 

Determination of the chemical composition of the gas in galax- 
ies versus cosmic time provides a very important tool for under- 
standing galaxy evolution, due to its important impact on fields 
such as stellar evolution and nucleosynthesis, gas enrichment 
processes, and the primary or secondary nature of the different 
species. 

Optical emission lines in galaxies have been widely used 
to estimate abundances in extragalactic Hn regions (e.g. Aller 
1942, Searle 1971,Pagel 1986, Shields 1990). The direct method 
of estimating metallicity in galaxies is known as the "T e method" 
(Pagel et al. 1992, Skillman & Kennicutt 1993), and consists of 
measuring the oxygen abundance from the electron temperature 
of the H ii region obtained using, the ratio of a high-excitation 
auroral line such as [O in] A4363 to the lower excitation [O m] 
AA 4959, 5007 lines. 

However, [O m] /14363 is too weak, not only in metal-rich, 
but even in metal-poor galaxies. In the absence of this auroral 
line, metallicities have to be estimated from strong line ratios, 
such as the R23 ratio inttoduced by Pagel et al. (1979). This 
is a commonly used and well-calibrated method of estimating 
metallicity abundances (see Alloin et al. 1979, Edmunds & Pagel 
1984, McCall et al. 1985, Dopita & Evans 1986, Pilyugin 2000, 
2001, Tremonti et al. 2004, Kewley & Dopita 2002, Liang et al. 
2006, hereafter L06). Nevertheless, this method has the disad- 
vantage of being double-valued as a function of 12+log(0/H). 
Besides R23, there are other strong-line methods useful for de- 
termining abundances of high-metallicity, star-forming galax- 
ies, for example the N2=[Nn] /16583/ Ha method (Denicolo et 
al. 2002, Pettini & Pagel 2004), the S 23 =([Sn] AA6111, 6731 



+[Sm] ,U9069,9532)/H/3 (Vflchez & Esteban 1996), among 
others. For a detailed discussion of the different methods see 
Bresolin (2006) and Kewley & Ellison (2008). 

Many studies of metallicity evolution exist as a function of 
cosmic time, although many of them refer to damped Lyman a 
systems at z > 2 (e.g. Pettini et al. 2002, Henry & Prochaska 
2007). For the evolution of the mass-metallicity relation of star- 
forming galaxies, Erb et al. (2006) find evolution at z > 2 and 
z ~ 3.5, respectively. Also, Brooks et al. (2007) and Finlator & 
Dave (2008), among others, derive cosmological models of the 
mass-metallicity relation, while at intermediate redshifts (1 < z 
< 2), there are several important studies of the evolution of the 
chemical composition of the gas, such as the ones by Lilly et al. 
(2003), Kobulnicky et al (2003), and Maier et al. (2006). 

Among the studies at z < 1, mainly based on small samples, 
Savaglio et al. (2005) have investigated the mass-metallicity re- 
lations using galaxies at 0.4 < z < 1, finding that the metallicity 
is lower at higher redshift, for the same stellar mass, by 0.15 
dex. Also, Maier at al. (2005), with a sample of 30 galaxies with 
0.47 < z < 0.92, find that one-third have metallicities lower than 
those of local galaxies with similar luminosities and star forma- 
tion rates. In contrast, Carollo & Lilly (2001) use emission-line 
ratios of 15 galaxies in a range of 0.5 < z < 1 to find that their 
metallicities appear to be remarkably similar to those of local 
galaxies selected with the same criteria. A similar result was 
found for the luminosity-metallicity relation by Lamareille et al. 
(2006), comparing star-forming galaxies at local and intermedi- 
ate redshift (0.2 < z < 1, split into 0.2 redshift bins). However, 
Buat et al. (2008) and Kobayashi et al. (2007) have derived a 
model of metallicity evolution as a function of z, which show a 
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progressive increase in metallicity with time, even at low red- 
shift. 

There is, then, a need to increase the lower redshift galaxy 
samples to ascertain whether or not at such low redshifts (i.e. 
small lookback times, of the order of 8.4 Gyr for z ~ 0.5, using a 
concordance cosmology) any evidence exists of metallicity evo- 
lution. The SDSS database provides an excellent opportunity for 
extending these studies up to z ~0.4 to explore a possible evo- 
lution of metallicity at low-redshift using larger samples, thus 
deriving more statistically significant results. 



2. Sample selection 

We used the SDSS-DR5 (Adelman-McCarthy et al. 2007) spec- 
tra from the STARLIGHT Database^, which were processed 
through the STARLIGHT spectral synthesis code developed by 
Cid Fernandes and colleagues (Cid Fernandes et al. 2005, 2007, 
Mateus et al. 2006, Asari et al. 2007). We obtained the emission 
lines fluxes measurements of our samples from the continuum 
subtracted spectra. For each emission line, STARLIGHT code 
returns the rest frame flux and its associated equivalent width 
(EW), linewidth, velocity displacement relative to the rest-frame 
wavelength and the S/N of the fit. In the case of Balmer lines, the 
underlying stellar absorption was corrected using synthetic spec- 
tra obtained by fitting an observed spectrum with a combination 
of /V* = 150 simple stellar population (SSPs) from the evolution- 
ary synthesis models of Bruzual & Chariot (2003). 

Since we are interested in studying the chemical evolution 
of emission-line galaxies at different redshifts, our sample was 
divided into three intervals: z\ =(0.1 -0.2), Z2=(0.2-0.3), and 
Z3=(0. 3-0.4). To avoid luminosity biases, we selected galax- 
ies in the same luminosity intervals, taking as a reference the 
magnitude completeness of the farthest interval Zi- Then, we se- 
lected sample galaxies with absolute Petrosian r magnitude in 
the range -24.8 < M r < -23.1 mag. As a consequence, we end 
with a sample spanning a narrow luminosity range, thus avoid- 
ing luminosity-metallicity evolutionary effects, as will be shown 
in Fig. 5. Absolute magnitudes were both, K and Galactic ex- 
tinction corrected by using the code provided by Blanton et al. 
(2003) and the maps of Schlegel et al. (1998), respectively, as 
provided by the STARLIGHT team. The Schlegel resulting sam- 
ple contains 25812 galaxies forz3, 15193 forz2 and 4225 forzi. 
It was not possible to extend the sample to redshifts lower than 
0. 1 due to the scarce number of spectra. 

From this sample we only consider galaxies whose spec- 
tra show in emission the Ha, U/3, [Nn] ,16583, [On] ,13727, 
[O m] ,14959, and [O ra] ,15007 lines, with a signal-to-noise ra- 
tio for [On] ,13727 higher than 3<x. A final selection for star- 
forming galaxies excluding AGNs was made using the criteria 
given by Kauffmann et al. (2003) in the [Om] i5007/H/3 vs. 
[Nn] /i6583/Ha' diagram, leaving a sample of 34 galaxies for 
Z3, 142 for Z2, and 40 for z\ ■ 

Since Balmer lines have already been corrected for underly- 
ing stellar absorption by the STARLIGHT code, it is only nec- 
essary to correct emission lines for dust extinction. Our extinc- 
tion correction was derived using the Cardelli extinction curve 
(Cardelli et al. 1989) and the Balmer decrements, assuming case 
B recombination for a density of 100 cirT 3 and a temperature of 
10 4 K, results in a Ha/H/? predicted ratio (unaffected by redden- 
ing or absorption) of 2.86 (Osterbrock 1989). 



1 http://www.starlight.ufsc.br 



3. Metallicity estimates 

We estimate metallicities using the R23 relation, introduced 
by Pagel et al. (1979), fl 2 3=([On] ,13727] + [Om] X14959, 
5007])/Hy6, adopting the calibration given by Tremonti et al. 
(2004), 12+log(0/H) = 9.185 - 0.313x - 0.264jc 2 - 0.321x 3 , 
where x - log R23 ■ 

However, this calibration is valid only for the upper branch 
of the double-valued R23 abundance relation, and additional line 
ratios, such as [N 11] ,16583/[0 n]/13727, are required to break this 
degeneracy. Since the upper and lower branches of the R23 cali- 
bration bifurcates at log([N n]/[0 n]) — 1.2 for the SDSS galax- 
ies (Kewley & Ellison 2008), which corresponds to a metallicity 
of 12+log(0/H) ~ 8.4, we select galaxies having 12+log(0/H) > 
8.4 and log([Nn]/[On]) > -1.2, corresponding to the upper R23 
branch. Applying this criterion, we end up with a final sample of 
28 galaxies for Zi, 140 for zi, and 39 for z\ , by discarding galax- 
ies of the lower branch we are not introducing a bias because 
96% of our galaxies lie in the upper branch. 

From these data, we derived the abundance-sensitive diag- 
nostic diagram [Nn] /l6583/[On] ,13727 vs. 12+log(0/H), rep- 
resented in Fig. 1. We selected this diagram due to its low scat- 
ter and physical information. The advantages of using [N n] and 
[0 11] lines is that they are not affected by underlying stellar pop- 
ulation absorption, and because this ratio is almost independent 
of the ionization parameter, since N + and + have similar ion- 
ization potentials. 
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Fig. 1. Calibration relation between 12+log(0/H) and log([Nn] 
/16583/[Oii] /13727) for the three samples in redshift, triangles 
represent galaxies of the sample at redshift z\, squares galaxies 
of Z2, and asterisks galaxies of zs- The solid line represents the 
best fit of the data using a linear fit, with the 2cr discrepancy 
in short dashed lines, and the dot-dash line for an order three 
polynomial fit given by L06 

This diagram has also been used by Kewley & Dopita 
(2002), Nagao et al. (2006), and L06 among other metallicity- 
sensitive emission-line ratios, like log([Nn] /16583/Ho'), 
log([Ora] /15007/H/3)/[Nn] ,16583/Ha), and log([Ora] ,14959, 
5007/H/?). Figure 1 shows our three samples in redshift, as well 
as the fits of L06 for their sample (0.04 < z < 0.25, and Petrosian 
r magnitude in the range 14.5 < m r < 17.77). 
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Fig. 2. Normalized metallicities histograms for our three sam- 
ples. The dashed line represents the sample in the redshift inter- 
val zi, the sample at zi is symbolized by a solid line, and the 
dot-dashed line corresponds to the sub-sample at redshift Z3 

To interpret our results, it is important to note that we are 
working with the integrated spectra, but that we are comparing 
with previous studies of nuclear spectra. As shown by Kewley 
et al. (2005) and L06, using a sample of Jansen et al. (2000), 
data points from nuclear spectra follow the SDSS galaxies nu- 
clear spectra very well, but data points from the integrated spec- 
tra show lower 12+log(0/H). The nuclear metallicities exceed 
the metallicities derived from integrated spectra by ~0. 13 dex 
on average. 

Based on this diagram, we claim that an evolution with red- 
shift is present in our sample . There is a clear decrement of 
~0.07 dex in 12+log(0/H), corresponding to the samples of zi 
and zi compared to the L06 polynomial fit, an effect resulting 
from the fiber size with respect to the size of the galaxies, as 
pointed out by Kewley et al. (2005). Both samples are also ~QA 
dex lower in [Nn] /i6583/[On] /i3727. However, the zi sample 
is ~0.1 dex lower in 12+log(0/H) with respect to the zi and zi 
samples. Samples corresponding to z\ and zi, seem to follow the 
same distribution in metallicity, but the Z3 sample shows a clear 
decrement in 12+log(0/H) of -0.1 dex and [Nn] ^6583/[On] 
/13727 of ~0.2 dex related to the other two samples. 

4. Metallicity evolution 

To validate the decrement observed in Fig. 1, we proceed to gen- 
erate a histogram of metallicities for our three redshifts samples 
as shown in Fig. 2. In the histograms we can appreciate a shift 
to lower metallicities as redshift increases, even for the z\ and Z3 
samples, in which we have only 39 and 28 galaxies, respectively. 

Although the derived distributions suggest a metallicity evo- 
lution, to investigate whether or not these could be an artifact due 
to the limited number of galaxies, we performed Monte Carlo 
simulations to determine the probability that the three distribu- 
tions represent the same sample. Taking as a reference the zi 
distribution, because it has the larger number of galaxies, we 
simulated this distribution by fitting a Gaussian, but using the 
number of galaxies of the other two samples. After a thousand 
simulations, we find an 8% and a 3.5% probability that this dis- 



tribution represents the same as that of z\, and Z3, respectively. 
This result supports the idea of an intrinsic evolution in metal- 
licity as observed at 23. 

5. Discussion 

To investigate whether the origin of the decrement in metallicity 
of sample zi is due to instrumental effects or to an inherent prop- 
erty of the sample of galaxies, it is necessary to explore three 
important effects: the origin of nitrogen, the effect of the 3 arc- 
sec diameter of the Sloan fibers, and the luminosity-metallicity 
relation. 

5.1. Nitrogen and oxygen abundances 

The observed decrement in metallicity at the higher redshift 
could be due to the primary or secondary origin of nitrogen as 
redshift increases. To determine the origin of nitrogen in our 
samples, we estimate the abundance for these galaxies, a de- 
tailed explanation of this method will be given in our next pa- 
per (in preparation). We do not appreciate a significant differ- 
ence between our redshift samples concerned with the primary 
or secondary origin of nitrogen. For the sample Z3, a 57.2% of the 
galaxies lie in the primary zone of the nitrogen, and the rest in 
the secondary part. Therefore, the decrement observed in metal- 
licity is not due to a bias in the primary or secondary origin of 
nitrogen in the samples. Our results agree with those obtained by 
Vila-Costas & Edmunds (1993), L06, Kennicutt et al. (2003). 

5.2. Effect of the Sloan fiber diameter 

With respect to the Sloan fibers, we expect that at higher redshift, 
the projected size of the Sloan fibers will cover a larger fraction 
of galaxy area than for nearby galaxies. As argued in Sect. 3, 
since integrated metallicities are lower than nuclear, this could 
explain the metallicity decrease that can be observed for the Z3 
sample with respect to the other ones. To quantify this effect for 
our sample, we estimated the percentage of angular size of each 
galaxy inside the Sloan fiber using the Petrosian total radius in r 
band in arcsec. 

We find that the 3 arcsec diameter of the Sloan fibers cover 
a maximum of 50% of the angular size of the galaxies for the 
100% of the zi sample, 98% of the zi sample, and 90% of the z 3 
sample. It is clear that the differences are not significant and that 
the observed decrement in metallicity should not be attributed to 
a differential projected size of Sloan fibers at each sample. 

5.3. Effect of the luminosity-metallicity relation 

We also discard the effects of a luminosity-metallicity relation, 
because the galaxies in our sample do not show a dependence 
of luminosity versus metallicity, as show in Fig. 5, since we are 
using the same range of 1.7 absolute magnitudes at every red- 
shift interval, which is too small to find significant differences 
in metallicity. Important variations in the luminosity-metallicity 
relation can only be seen when spanning ranges of ~5 magni- 
tudes. 



6. Conclusions 

Although similar studies of metallicity using SDSS exist, they 
are either restricted to a redshift ~0.1 (e.g. Tremonti et al. 2004) 
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Fig. 3. Luminosity-Metalicity relation for our sample of galax- 
ies, symbols follow the same code used in Fig. 1. 



or do not separate their samples as a function of redshift (e.g. 
L06), thus making it impossible to detect metallicity evolution. 
Finally, studies at high redshift are statistically limited. 

After exploring possible biases on the sample, the results ob- 
tained in the present work suggest prima facie evidence of an 
intrinsic metallicity evolution in the local Universe, showing a 
decrement of ~0.1 dex in 12+log(0/H) at a redshift interval 0.3 
< z < 0.4, a result consistent with the models of metallicity evo- 
lution in Buat et al. (2008) and Kobayashi et al. (2007). 
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